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Abstract. The totally asymmetric simple exclusion process (TASEP) is a well 
studied example of far-from-equilibrium dynamics. Here, we consider a TASEP with 
open boundaries but impose a global constraint on the total number of particles. In 
other words, the boundary reservoirs and the system must share a finite supply of 
particles. Using simulations and analytic arguments, we obtain the average particle 
density and current of the system, as a function of the boundary rates and the total 
number of particles. Our findings are relevant to biological transport problems if the 
availability of molecular motors becomes a rate-limiting factor. 
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1. Introduction 

Over the last two decades, the study of driven diffusive systems [H [21 [3l H] has formed 
a significant branch in the general pursuit to understand non-equilibrium statistical 
mechanics. The underlying dynamics violates detailed balance, so that the long-time 
behavior is governed by a true non- equilibrium steady state (NESS). Highly unexpected 
and non-trivial properties are manifested, even in one-dimensional systems with purely 
local dynamics. A paradigmatic model in the last class is the totally asymmetric simple 
exclusion process (TASEP) [51 El [3, El EJ Hj. Characterized by open boundaries and 
particle transport, it displays the key physical signatures of a system driven far from 
equilibrium. In particular, its NESS develops a complex phase diagram controlled by 
the interactions of the system with its environment. While an understanding of these 
properties is of fundamental theoretical interest, the TASEP and its generalizations have 
also acquired fame as models for practical problems in traffic flow [TQl [TT] and biological 
transport [121 [131 HS1 ESI dZl HBl HSl 1201 EI]. 

In this paper, we study a simple variant of the basic TASEP, by introducing a 
global constraint on the available number of particles. Deferring the motivations for 
such a model to the next paragraphs, let us briefly summarize the essentials. Each site 
of a one-dimensional lattice is either empty or occupied by a single particle. Following 
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random sequential dynamics, the particles enter the lattice at one end (e.g., the "left" 
edge) with a given rate a, hop to the right with rate 7 (scaled to unity), subject to 
an excluded volume constraint, and exit at the far end with a rate (3. In the standard 
version of the model, both a and (3 are constant rates, independent of the number of 
particles on the lattice. Thus, we may regard the lattice as being coupled to a reservoir 
- or a pool - with an arbitrarily large number of particles. Here, we consider a TASEP 
with a finite supply of particles and report our numerical and analytical findings. A 
fixed number of particles, N to t, is shared between the lattice and the reservoir, in such 
a way that the entrance rate, a, depends on the number in the pool: N p . As a result, if 
more particles are found on the lattice, the number of available particles in the pool is 
lowered, leading to a corresponding decrease of a. Since the number of particles on the 
TASEP lattice, N, feeds back into a through the constraint N + N p = constant = N to t, 
we will refer to our model as "a constrained TASEP." For simplicity, we assume the 
reservoir to be so large that particles exiting the lattice are unaffected by N p , leaving 
j3 unchanged. Our focus here is to explore the consequences of limited resources: How 
will the phase diagram and the properties of the phases be affected as the total number 
of available particles is reduced? 

This behavior mimics the limited availability of resources required for a given 
physical or biological process. For example, in protein synthesis, ribosomes bind near the 
start sequence of a messenger RNA, in order to translate the genetic information encoded 
on the RNA into the associated protein. Ribosomes are large molecular motors which are 
assembled out of several basic units, and numerous ribosomes can be bound to the same, 
or other, mRNAs, so that multiple proteins are translated in parallel. When a protein 
is completed, the ribosome is disassembled and recycled into the cytoplasm. Under 
conditions of rapid cell growth, ribosomes or their constituents can find themselves in 
short supply, so that a self-limitation of translation, mimicked via a modification of 
a, can occur. In a previous study [H], a different aspect of ribosome recyling was 
considered. This work focused on the enhancement of the ribosome concentration at 
the initiation site, as a result of diffusion of the ribosome subunits from the termination 
site (of the same mRNA). Due to the spatial proximity of initiation and termination 
sites, a is effectively increased. By contrast, our investigation of constrained resources 
leads to an effective reduction of a. 

In the context of traffic models, our problem corresponds to a generalization of the 
parking garage problem [22]. Here, one special site (the "parking garage", or reservoir) 
is introduced into a TASEP on a ring (lattice with periodic boundary conditions) and, 
for this site only, the occupancy is unlimited. Particles ("cars") jump into the garage 
with unit rate (corresponding to (3 = 1), irrespective of its occupancy. Particles exit 
the garage with rate a, provided the site following the garage is empty. As we will see 
presently, this is a special case of our more general model. 

We begin with a description of our model and its observables of interest. Next, we 
outline results from simulations and compare them to a simple theory which builds on 
exact analytic results for the standard TASEP. We find excellent agreement for nearly 
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all (a,/3) and conclude with some comments and open questions. 
2. The model 

The standard TASEP is defined on a one-dimensional lattice of length L, with sites 
labeled by i, i — 0, 1, L — 1. Each site can be empty or occupied by a single particle, 
reflected through a set of occupation numbers {n^} which take the values or 1. The 
boundaries of the lattice are open. Particles hop from a reservoir onto the lattice with a 
rate a. Once on the lattice, particles will hop to the nearest-neighbor site on the right, 
provided it is empty. Once a particle reaches the end of the lattice, it hops back into 
the reservoir with rate (3. The dynamics is random sequential, leading to fluctuations 
in the local occupations as well as in the total number of particles, defined via 



The overall density is given by p = N/L. Ensemble averages are denoted by (...). 

Three distinct phases can be identified and are displayed in an (a, (3) phase diagram 
(see Fig. [1]): a low density (LD), a high density (HD), and a maximum current 
(MC) phase. The phases are distinguished, respectively, by their average densities: 
p = a = p_, p = 1 — (3 = p + , and puc — 1/2- The corresponding stationary 
currents are given by J_ = a(l — a), J + = (3(1 — (3), and Jmc = 1/4. In all of 
these expressions, finite-size corrections have been neglected. Because of particle-hole 
symmetry, all aspects of the HD and LD phases are related. The phase boundaries 
between the HD and MC phases and between the LD and MC phases mark continuous 
transitions. The line separating the HD and LD phases, a = (3 < 1/2, is a coexistence 
line. Here, the system consists of a region with low density followed by one with 
high density (p+), connected by a microscopically sharp "shock." This shock diffuses 
freely between the ends of the system, so that the average density profile is linear. Often 
this is referred to as the "shock phase" (SP). 

The system of interest here differs from the standard TASEP in one important 
way: The number of particles in the reservoir, N p , is finite, and we choose the on-rate 
a to depend on N p as follows. So as to distinguish this varying rate from the a of the 
ordinary TASEP, we denote the former as the effective on-rate a e ff and write 



where f(x) is a function that satisfies three conditions: (i) /(0) = 0, (ii) /(oo) = 1, 
and (iii) f(x) is monotonically increasing. The first of these conditions is self-evident; 
the second simply connects our model to the standard one with unlimited resources; 
the last is just common sense. In this notation, the "parking garage" problem [22] is 
characterized by f(N p ) = 0(N p ) where 9 is the Heaviside step function. Here, we wish 
to model, say, a cell with finite number of ribosomes, so that it is natural to assume 
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a eff ( N p) = af{N p ) , 
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Figure 1. Phase diagram of the standard TASEP. 



that a smaller number of particles in the pool results in a lower on-rate, etc., and that / 
would be a smoother function. Specifically, for the simulations shown below we choose 

f(N p ) =tanh(N p /N*) (1) 

where N* provides the scale for crossover to saturation. The detailed choice of iV* is 
unimportant, but it is convenient to make it extensive in L. For easier comparison, we 
choose iV* to be the average (stationary) density of the unconstrained TASEP defined 
by the parameters (a, (5): For example, if a > 1/2, (3 < 1/2 (HD phase, ordinarily), we 
set iV* to Lp + = L(l — (3). For the SP, we arbitrarily choose N* = L/2. To summarize, 
the control parameters of our model are the same as the ordinary TASEP (a, /3; L) plus 
the total number of particles, N tot (i.e., N + N p ). 

To characterize our system in steady state, we measure the average density 
p = (N)/L, the average (local) density profile pi = (rii)/L, and the average current 
J. For the last quantity, we record the total number of particles entering and leaving 
the chain over the run (thereby improving statistics) and divide the average of these by 
the length of the run. 

In our Monte Carlo simulations, we pick randomly from the N particles on the 
lattice and one additional virtual particle. If a lattice particle is selected, we attempt to 
update its position as in ordinary TASEP; if the virtual particle is selected, we attempt 
to place a new particle on the first site of the chain. N+l attempts constitute one Monte 
Carlo Step (MCS). Initially, the chain is empty, i.e., N = 0, N p = N tot . Typically, the 
first 10 6 (10 7 for SP) MCS are discarded, to ensure that the system has reached steady 
state, before data are taken. Then, configurational observables are recorded every 100 
MCS and averaged until the run terminates. The typical length of each run is 10 7 (10 8 
for SP) MCS. Longer runs are required to obtain low noise data for density profiles and 
currents. Statistical errors were estimated through visual inspections of time series and 
through multiple runs, to ensure the reproducibility of the data. We simulated system 
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sizes in the range of L = 250 to L = 8000, with most data taken for L = 1000. To 
obtain the phase diagram, we simulated more than a dozen (a, (3) pairs. Here, we report 
results for four points: (0.25,0.75), (0.75,0.25), (0.75, 0.75), and (0.25, 0.25). Since these 
correspond, respectively, to LD, HD, MC, and SP in a standard TASEP, we will use 
these letters to refer to the four cases studied. 

3. Results from simulations 

In this section, we summarize our simulation results. Since only the on-rate depends 
on the feedback, the behavior of our system is no longer governed by particle-hole 
symmetry. As a result, the simplest case is the first in the above, i.e., "LD". Both p 
and J increase with N tot in an expected fashion. The next case, MC, is already slightly 
more interesting: Both p and J increase until the effective on-rate reaches 1/2, after 
which both become constant. The end result is a pronounced kink in p. The HD case 
provides an even more interesting scenario: The system properties range through three 
regions as N to t is increased. Finally, the SP provides the most unexpected behavior. In 
the next section, we turn to a simple theoretical description of our findings. 

The LD case. Here, we set a < 1/2, (3 > a, and vary the total number of particles, 
N tot , in the system. With / given by Eqn ([1]), our simulation results for L = 1000 
and (a,/3) = (0.25,0.75) are shown in Fig. [2j As expected, we see that, for large 
Ntot, the system density p = (N) /L approaches the value for the standard TASEP for 
Ntot ^ N* (here, 0.75L), namely, p = a (which is 0.25 in our case). More interestingly, 
for N to t < N*, we find a reduced density, p < a, a signal that the system is responding 
to the limited availability of particles from the reservoir. In the limit N tot — ^ 0, p 
naturally vanishes, with a predictable slope (see below). The current J in this "phase" 
also decreases monotonically with decreasing N to t, from its asymptotic limit (3(1 — (3) 
(given by 0.1875). 
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Figure 2. Average density p and current J vs Ntot/L for an LD case (L = 1000, 
a = 0.25, and f3 = 0.75). Open circles and diamonds are from the analysis in Section 
4. 



Far- from- equilibrium transport with constrained resources 



6 



The MC case. This domain of the phase diagram is characterized by a > 1/2 and 
(3 > 1/2 in the standard TASEP. Fig. [3] shows the averages p and J, as N tot varies, 
for (a,/?) = (0.75,0.75). The former displays a single kink, just below N tot /L ~ 1.0, 
accompanied by a rather smooth crossover in the current. As we will see below, this 
kink is associated with a crossover from an MC-like behavior to an LD-like behavior, as 
N tot drops below the value required to sustain a density p MC = 1/2 on the chain. 




Figure 3. Average density p and current J vs N to t/L for an MC case {L = 1000, 
a = 0.75, and (3 = 0.75). Open circles and diamonds are from the analysis in Section 
4. 




Figure 4. Average density p and current J vs N to t/L for an HD case (L = 1000, 
a = 0.75, and (3 = 0.25). Open circles and diamonds are from the analysis in Section 
4. The predictions from domain wall theory (not shown) are similiar. 

The HD case. Here, we set (3 < 1/2 and a > (3. Fig. H] shows our results 
for (a,/3) = (0.75,0.25). As expected, for sufficiently high Ntot, the system settles 
into the density and current associated with the HD phase of the standard TASEP: 
p = 0.749 ~ p + , and J = 0.184 ~ (3 (1 — (3). As N tot is reduced, to the extent that 
the chain is prevented from sustaining a density of p + , the naive expectation is that a 
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crossover to an LD-like behavior appears. Remarkably, however, Fig. H] shows not just 
two, but three distinct regimes, separated by two pronounced "kinks" at N tot /L ~ 0.5 
and 1.0. As Fig. [5] illustrates, these kinks become sharper with increasing L. Looking 
for potentially critical behavior, we measured the variance of (N) near the kinks and 
compared it to its values deep within each regime. To our surprise, at least for the 
L = 1000 system, the variance was not noticeably larger at the kinks than at other 
locations. Simulations at larger values of L would be required to settle this issue with 
more certainty. 

The current J is also displayed in Fig. HI For low values of N tot /L, J increases 
in an expected way. Around N to t = 0.5L, it shows a distinct kink and then becomes 
essentially constant, with a value of 0.184. This is close to 0.1875, the value of J for 
infinite L. Since the current appears to have already reached the "asymptotic" value 
here, it is not too surprising that there is no sign of a "second kink" (near N to t = L, 
where p has another one). Below, we offer a theoretical description for these findings. 
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Figure 5. Finite-size effects on the kink in p, near Ntot/L — 1.0, for the HD case in 
Fig. 4: L = 250, 1000, and 4000. 

The SP case. Defined in the standard TASEP by a = (3 < 1/2, this line of first- 
order phase transitions separates the HD and LD phases. Though the average density 
is 1/2, its fluctuations are large. In particular, since a shock (which separates a region 
with density p + from one with p_) diffuses freely between the ends of the system, the 
probability P(p) of finding the system with a density p approaches a flat distribution: 
P(p) = 9(p — p-)6(p + — p)f (p+ — p_), in the limit L — ► oo. Of course, finite-size effects 
will smooth out the steps at p + and p_. Imposing finite resources, we may expect 
that p increases simply with N tot as in the LD case, perhaps rising more slowly to the 
asymptotic value of 1/2. Instead, its properties are much more subtle. As illustrated in 
Fig. M (for a = P = 0.25), p increases with N tot /L, slows, speeds up, and slows down 
again. The behavior in the low N tot regime is intuitively understandable, since there 
are too few particles to sustain any phase other than LD. The large Ntot regime is also 
expected. In the crossover regime (L < N tot < 3L), an intriguing "shoulder" develops. 
Though this behavior is accounted for by the domain wall theory [231124"] (see below), we 
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have found no simple and intuitive way to understand it. Perhaps the presence of large 
fluctuations, as the system gets closer to the SP, smooths out the two kinks observed in 
the HD case. To gain more insight, we investigated P(p). Fig. [7] shows particle density 
histograms (gathered from 50K data points and plotted with arbitrary normalization 
here), for several choices of N tot /L. For low N tot /L (e.g., 0.5 in the figure), our system is 
dominated by the lack of resources, so that this distribution is essentially Gaussian, quite 
similar to those measured deep within the LD phase. As N tot /L increases, the system 
enters the crossover regime and P{p) becomes quite asymmetric, as a detailed fitting of 
the N tot /L = 1.0 case shows. At the other extreme, this distribution is indeed relatively 
flat. However, the true asymptotic regime is reached only for very large N tot /L. Thus, 
the slope for even the N tot /L = 3.0 case (cf. Fig. [7J is discernably non-zero. Throughout 
the intermediate regime, the structure of P is more interesting. Each histogram has the 
cross-section of a "lean-to" (a shack with a slanted roof). A more detailed inspection 
shows that the "roof" is predominantly a (slowly) decaying exponential 0. This picture 
is consistent with the one from domain wall theory [23j El], as will be discussed below. 
It is this slow crossover - from well-defined Gaussians to a completely flat distribution 
- which accounts for the interesting structure in the rise of p with N to t/L. 

Finally, the behavior of the current is also sensitive to the slow crossover, although 
it displays no interesting structures like those in p. Comparing Fig. [6] with the earlier 
ones, we see that it saturates at values of Ntot/L which are 2 to 3 times higher than in 
the other cases. 




Figure 6. Average density p and current J vs Ntot/L for an SP case (L = 1000, 
a = 0.25, and (3 = 0.25). Open circles are from the domain wall theory in Section 4. 
Open diamonds are the current computed from Eqn (|2j). 



| Such distributions are also found in unconstrained TASEP's with a nearly identical to (3. 
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Figure 7. Histograms of the density p for L = 1000, a = 0.25, and j3 = 0.25. (For ease 
of comparisons, the histograms with solid lines only are scaled down by 2.) Legend 
shows various values of N to t/L. 



4. Analytic predictions 

In the steady state of a standard TASEP, the overall density of particles is controlled 
only by the entrance and exit rates. It is natural to expect that these relations would 
also hold in our case, except that the stationary density needs to be determined from 
the (variable) entrance rate in a self-consistent way. In the following, we pursue this 
approach and compare it to the simulation results reported in the preceding section. 

The LD case. Again, we begin with a choice of (a, (3) such that the system is in an 
LD phase, with p < 1/2. Invoking the exact expression for the density of the standard 
TASEP, namely, p = a, we write the density of particles in the resource-limited 
p = a e ff = af(N p ). With N p = N tot — Lp, this becomes an implicit equation for p: 

P = af(N tat - Lp) (2) 

which can be solved for given a and N tat . Specifically, since < f(N tot — Lp) < 1, the 
solution will lie between and a. The upper bound will be reached if N tot /L —* oo so 
that the right hand side of Eqn (j2j) remains essentially constant as a function of p. For 
any finite value of Ntot/L, the right hand side is a monotonically decreasing function of 
p so that (i) Eqn (j2J) always has a solution, and (ii) the solution decreases monotonically 
with Ntot/L. The resulting density, and the associated current, J = a e ff (1 — « e //), are 
shown as open circles and diamonds, respectively, in Fig. [2J for the choice of / given in 
Eqn (Til). We see, in particular, that for small values of N to t/L, we may use tanhx ~ x 
whence p ~ N tot /(2L). Clearly, the agreement with the simulation data is impressive. 

The MC case. After these considerations, the analysis of the MC phase is quite 
straightforward. For large N to t, the density saturates at p = 1/2, with J = 1/4. For 
small N to t, the system is LD-like and follows Eqn (j2J). The transition occurs when a e ff 
is forced to drop below 1/2, i.e., at the intersection of Eqn (j2J) and p = 1/2. Here, N to t 
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assumes the critical value of N totiC , given by 

AW = Lp + r 1 (±) . (3) 



2a, 

The current follows the behavior of the density. 

The HD case. Next, we turn to the high-density phase, with (3 < 1/2 and a > (3. 
If the number of available particles is very large, we may again invoke the exact relation 
for the density of the standard TASEP, namely 

p = p+ = \-p. (4) 

This reflects the constant density observed in the simulations. Clearly, this cannot 
be sustained as the number of available particles is reduced, since the effective a 
of the system is also reduced. A transition should be expected when it reaches (3. 
= a e ff = af(N p ). Translating (3 = a e ff into an equation for the system density p, we 
arrive at 

P = a eff = af(N tot - Lp) . (5) 

Solving for p, and setting it to the one from Eqn (jlj), we find the critical value of the 
total particle number for this transition: N to t.2 = L (1 — (3) + f^ 1 ((3/ a). (The subscript, 
2, helps to distinguish it from the lower transition point.) For the specific / used in 
simulations, this results in 



N tot . 2 = L{\-(3) 



1 + tanh" 1 - 

a 



(6) 



As we see, Ntot.2/ L = 0.75 [1 + tanh" 1 (1/3)] ~ 1.01 is in excellent agreement with the 
"second kink" shown in Fig. HI 

At the opposite extreme, the resources are very low, so that a e ff is significantly 
less than (3, and the system is pushed into LD-like behavior. For such low N to t, the 
density follows Eqn (TS]). In this regime, / is well approximated by a linear function, so 
that p ~ a(N tot - Lp)/N* = a(N tot /L - p)/ (1 - p), i.e., 

P - 1 Ntat/L . 

1 — jj + a 

Of course, Eqn (J2J) is also not "sustainable," in the sense that p — > a for large N tot 
instead of the correct limit: p — > 1 — (3. The critical value of N tot again occurs when 
a e ff reaches /3, i.e., when the density given by Eqn (|2J) reaches (3. The result is 

N totA = LP + L (1 - (3) tanh" 1 ( . (7) 



Again, for the case shown in Fig. HI there is excellent agreement between the data and 
this prediction: N to t.i/L ~ 0.51. 

Between N tot .i and N tot _2, Eqn (jSJ) provides a linear dependence of p on N tot /L 

p = N^L - r 1 (PI a) /L , (8) 

with unit slope. Again, with no fitting parameters, this theoretical prediction agrees 
remarkably well with the observations. In this regime, the reservoir occupation (and 
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so cteff) remains constant, while the change in N tot is balanced by the change in total 
lattice population: Lp. Since a e ff = [3 (on the average), the current remains at the 
constant value J = /3 (1 — (3) throughout. 

While this "three-piece" approach is both intuitively appealing and surprisingly 
successful in reproducing simulation data, it begs the question: Is there a unified theory 
that is just as successful with predictions? The answer is affirmative, namely, the 
domain wall theory [231 |21| ■ Since this theory is absolutely indispensable for providing 
good predictions for the next case, we will discuss its details in the next subsection. 

The SP case. Here, Monte Carlo simulations presented us with unexpected 
phenomena: the "shoulder" in Fig. [6] for p. The simple approach taken above succeeds 
only in predicting the behavior for the extreme values of N tot /L. Dominated by the 
limited resources, the solution of p = a e ff fits the data very well for N tot /L < 1. At 
the other extreme (N to t/L ^> 1), the system is necessarily in the SP, where p = 1/2 on 
the average. In the crossover regime, the understanding of the structure in p(N tot /L) 
requires a more sophisticated theory [23, 24J. Based on a biased random walk of the 
shock (i.e., domain wall), Santen and Appert [21] derived the average occupation for a 
finite TASEP: 



where 5 = p + — p_ = (1 — (3) — a, and 



is the ratio of the diffusion constants to the right/left. To apply to our SP case, we 
simply set a to be a e ff = (3f{N tot — Lp) in this equation. That sets up a relationship 
between p and N tot /L which can be computed numerically. The agreement between 
this theory and data is quite respectable (Fig. [6]), though not as spectacular as in the 
other cases. Not surprisingly, most of the disagreement occurs in the crossover region. 
We believe that the main difficulty lies with the large fluctuations of the domain wall, 
feeding back in non-trivial ways to the on-rate. 

For another perspective on the successes and limitations of this theory, we can 
compare P(p) above with the predicted Pi, the probability of finding the domain wall to 
be at site i. The latter is, for the steady state, a pure exponential [23,121]: Pi oc exp (i/£). 
Since our focus here is the dependence on p, Pi can be translated into P(p), via the 
relation p = p + (1 — i/L) + p_ {i/L). We find that the predicted P(p) is proportional 
to exp (— p/er), and nonzero only in the interval [p_,p + ]. Our simulations show several 
non-trivial deviations from this theoretical prediction, however, (i) The slopes in the 
exponent are systematically steeper than predicted, (ii) Quadratic terms in In P(p) are 
far from negligible, (iii) The tails of P(p) extend substantially outside the range [p_, p + ]. 
In other words, there is room for improvements to the domain wall theory. 

Another intriguing phenomenon, not shown here explicitly, is the complex 
dependence of this crossover on L. In particular, both simulations and theory show 




(9) 
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that the "shoulder" completely disappears for small L. For example, in an L = 100 run, 
p{Ntot/L) rises linearly for N tot /L < 1.5 and simply "slows down" into the asymptotic 
value of 1/2 for N to t/L > 2.5. For large L, in contrast, there exists a subtle competition 
between the factors {L + 1) and r L+1 in Eqn (Q. The result is that the crossover 
region gets longer with larger L, producing challenges at the simulations front as 
well. Of course, these difficulties are hardly surprising, since the complexities here 
are undoubtedly due to large fluctuations associated with the infrared scales. Further 
studies into this issue, especially more detailed investigations of the lean-to's in P(p), 
should provide insight for a better understanding of such fluctuations. 

To conclude this section, we briefly note that we also simulated an alternate form of 
the on-rate, namely, f(N p ) = 1/ [l + (N*/N p ) ]. We will not present any details here, 
but the agreement of simulation data and analytic results for the three regimes (HD, 
LD, and MC) is equally convincing. However, the details of the crossover regime in SP 
are likely to be quite sensitive to functional form of /. 

5. Conclusions 

In this paper, we studied a TASEP with open boundary conditions. In contrast to 
the standard model, where the particles are supplied by an infinite reservoir, the total 
number of available particles in our model, N tot , is fixed. Mimicking limited resources, 
the available particles are shared between the reservoir and the system itself: If more 
particles bind to the chain, the reservoir is depleted, and vice versa. We account for 
this situation through an effective on-rate, a e //, which depends on N p , the number 
of particles in the pool, and hence, on N, the number on the chain. Specifically, we 
mainly used Eqn (Op), a e ff = atanh (N p /N*), here. As the number of available particles 
is reduced from infinity, the system crosses over from the asymptotic behavior of the 
standard TASEP to different behaviors dominated by limited resources. 

If the corresponding standard model (i.e., N to t °o limit) lies in the LD phase, 
there are no surprises in either p, the average particle density on the chain, or J, the 
average current: Both are proportional to N tot for N tot <C L, crossing gently over to the 
asymptotic values much like the tanh function. For the MC phase, there is already a 
sublety: both p and J increase monotonically and then, at a characteristic a e ff (= 1/2) 
become constant, p displays a notable kink there. For the other two cases (HD and SP), 
p exhibits an additional, intermediate regime (though the behavior of J shows little hint 
of this regime). In the HD case, p is linear in Ntot/L with unit slope; only the intercept 
depends on a/j3. Here, N p remains constant, as any increase in N tot /L is "absorbed" 
by the chain (through a shift in the domain wall). This behavior may be regarded 
in the same light as an equilibrium system with phase co-existence. Considering, e.g., 
the pressure vs. volume isotherm for a liquid-gas system below criticality, the pressure 
remains constant over some range, as any increase in volume is absorbed by a shift in 
the liquid-to-gas ratio. 

The most interesting and challenging case is the SP. In the intermediate regime, 
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a "shoulder" appears in p(N tot /L). In addition, this structure displays a non-trivial 
dependence on L. Although the predictions of domain wall theory agree reasonably well 
with MC data, there are small (< 6%) systematic discrepancies. We also investigated 
histograms of particle densities, which cross over from sharply peaked Gaussians for 
N to t <^ L to a flat distribution (between p_ = j3 and p + — 1 — (3) for N tot ^> L. In the 
intermediate regime, this distribution displays a much richer behavior than a truncated 
exponential - the result of domain wall theory. 

Several questions remain open for further study. First, are the transitions between 
the three regimes in HD associated with true thermodynamic singularities? We find that 
the changes in slope accompanying these "crossovers" become sharper with increasing 
chain length L. While this speaks for a true transition, the absence of large fluctuations 
is somewhat surprising. Clearly, a more careful analysis is required to explore these 
questions more fully. 

Second, though domain wall theory is reasonably successful for p(N tot /L) in the 
SP case, we would like to find improvements of this approach, so as to account for 
the crossover regime better. Arriving at a good understanding of P(p) and the large 
fluctuations would be desirable. Progress along these lines may also help us to develop 
an intuitive picture for the complex L dependence. 

Finally, we are mindful of one of the motivations of this study, namely, finite 
resources in biological systems, e.g., ribosomes (and aa-tRNA's) for modeling protein 
synthesis in a cell [121 EH (HI EEl CEO EH] ■ F° r that application, the above TASEP needs to 
be generalized, to include exclusion at a distance (particles covering more than one site) 
and inhomogeneous hopping rates. How these systems are affected by finite resources 
will pose many interesting new challenges, especially on the theoretical frontier. 
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